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Abstract 

The interaction of bodies in a fluid, mediated by hydrodynamic fluctuations 
and proposed by Dzyaloshinskii, Lifshitz, and Pitaevskii, is calculated ex- 
actly for parallel infinite planes and is shown to be attractive. The second 
mechanism of fluctuation interaction, proposed in the present work, is due 
to fluctuations of linear and angular velocities of bodies in a hydrodynamic 
medium and leads to a repulsion. The both mechanisms provide an inter- 
action energy of two particles in a fluid of the order of temperature for the 
inter-particle distance of the micron scale, where the interaction mediated by 
electromagnetic fluctuations is small. When two particles approaches a wall, 
placed inside a fluid, the first (attractive) interaction is violated a little, but 
the second (repulsive) weakens. This behavior correlates with the experimen- 
tally observed attraction of particles which appears when they approach a 
wall. 
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I. INTRODUCTION 



The interaction of colloidal particles is a matter of challenge and controversy despite of 
a long history of the problem An attraction of micron size particles separated by a 

distance of the micron scale has been established experimentally 0-0. An existence of the 
observable attraction at such big distance sets a real puzzle. Colloidal particles accept in a 
fluid a surface charge screened by counterions at Debye's length and an important question 
is that can those particles with like charges attract each other due to solely electrostatic 
mean-field interaction. This type of attraction was predicted in Ref. 0. The results of 
Ref. 1^ are not applicable to dielectric particles but only to those containing inside the 
electrolyte identical to outside one. The correct calculation of the interaction of that kind 
shows only repulsion in contrast to Ref. 0. The electrostatic attraction between like- 
charge particles based on the mean- field approach was also predicted in Ref. but that 
conclusion was incorrect, as shown in the work |Tl|. An attraction of like charged particles 
due to a mean-field mechanism seems to be extremely unlikely and the work [|ll|] provides 
strong arguments in support of this statement. The conventional van der Waals attraction, 
mediated by the high-frequency (visible light) electromagnetic fluctuations, is very small at 
the micron distance consisting less than 10~^T at room temperature |]T2[ (see also Ref. [0). 
Calculations of the effect of fluctuations of counterions |T^-|I^ is actually an account of some 
sort of the van der Waals interaction mediated by low frequency (generic with the plasma 
frequency of the electrolyte) electromagnetic fluctuations. It is remarkably, the attraction 
found in the above works is sufficiently short range (compared with the micron length) and 
it is much bigger than temperature what enables to compete with the Coulomb repulsion on 
such short distance. This contrasts to the experimental situation , where the attraction 
appears on micron distance when the Coulomb repulsion becomes of the order of temperature. 
As one can conclude, neither an electrostatic mean-field interaction nor the van der Waals- 
type attraction, mediated by high-frequency [|12| and low-frequency |T^-|1^] electromagnetic 
fiuctuations, can explain the observed attraction |^-|^ of two colloidal particles separated 
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by the micron scale distance. 

The first goal of this work is to calculate exactly the interaction of two parallel infinite 
planes mediated by hydrodynamic fluctuations and to draw a conclusion about the anal- 
ogous interaction of two particles. This type of fluctuation forces was discussed firstly by 
Dzyaloshinskii, Lifshitz, and Pitaevskii [TB|. It turns out, the interaction mediated by hy- 
drodynamic fluctuations is attractive and not small at the micron distance and is relevant 
for explanation of the observed attraction between colloidal particles The second 

goal of this work is to propose a new type of interaction of particles mediated by thermal 
fluctuations of linear and angular velocities (body variables) of particles in a hydrodynamic 
medium. This interaction is repulsive in contrast to one mediated solely by fluid fluctuations 
and in combination with it can explain the experimentally observed interaction of colloidal 
particles. 



II. HOW TO ACCOUNT THE HYDRODYNAMIC VAN DE WAALS 

INTERACTION 

Suppose two bodies are placed inside a hydrodynamic medium, they are totally fixed 
in space, and serves only as obstacles for a fluid motion due to stick boundary conditions 
at their surfaces. There is no any macroscopic motion in the system and the only motion 
is caused by thermal fluctuations of the fluid velocity v{f,t). In this case the free energy 
of thermal fluctuations of the fluid F{R) depends on the distance R between bodies. The 
function 

U,dw{R) = F{R) - F{^) (1) 

is an interaction mediated by fluid fluctuations. Analogously to conventional van der Waals 
interaction mediated by electromagnetic fluctuations, the potential (|T]) can be called the 
hydrodynamic van der Waals interaction. To find the free energy of thermal fluctuations of 
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the fluid one can start with the hnearized Navier-Stokes equation |T9 



= - + r/V^^T + (C + I) Vt; (2) 
The solution of Eq. |^ can be written in the form 

V = curl a + (3) 
where the vector a satisfies the equation 

da ^r,^ 

PQ^ = ^V^a (4) 
The scalar is determined by the equation 

,^._,^vV+fc + ^?VW) (5) 



dt ' y 3 ^ 

In Eq. ^ the schematic expression of Vp through the sound velocity c = ^Jdjyjdp is used. 
The density obeys the continuity equation dp/dt = — pV^0 what, together with Eq. 
provides the linear equation for V^0. Eq. ^ sets the typical frequency of fluctuations rjk'^/ p, 
where A; is a wave vector. At this frequency the term c^V^p ~ in Eq. |^ dominates for 
sufficiently small wave vectors k < I/Rq, where 

Ro = - (6) 

pc 

is some characteristic distance. The domination of that term requires the condition 

V'^ = (7) 

From the physical point of view it means that on distances bigger than Rq the fluid fluctua- 
tions correspond to an incompressible medium and fluctuations of the longitudinal velocity 
can be neglected. The vector a can be expanded by means of eigenfunctions tjji{f) 

i 

of the boundary problem formulated in the following way 

— \/'^a = pd, (curl a + V0) |b = (9) 



where B denotes the total boundary of all participating bodies and satisfies Eq. 0. Eq. ^ 
now becomes of the form 

Here Sj are some independent displacements of the fiuid chosen according to the above 
diagonalization procedure. Now one should write down the free energy associated with the 
set (|10]) of independent harmonic oscillators, what is done in the next Section. 

III. FREE ENERGY OF THE DAMPED OSCILLATOR 

Let a harmonic oscillator obeys the classic dynamic equation 

5^ + 75^ + tf. = (11) 

Due to thermal and quantum oscillations the free energy of this oscillator is not zero. In 
the under-damped limit 7 = the expression of the free energy T In (2 sinh(fin/2T)) is 
well known. For non-zero 7 a calculation of the free energy constitutes a typical problem 
of the field called dissipative quantum mechanics. The initial ideas of this field, steamed 
to Feynman [^, have been essentially developed by Caldeira and Leggett and other 



researchers, see the book ||22| and references therein. One of the way to calculate the free 
energy associated with the damped harmonic oscillator (0) is to consider a harmonic elastic 
string attached to the free particle in the potential mVt^s^ /2. In the classic limit this provides 
the equation of motion ([TTD with the finite 7 due to a radiation friction. Then the total 
harmonic system "particle and string" can be diagonalized to determined eigenfrequencies 
uji which differ from eigenfrequencies uo^ of the string without a particle. The free energy of 
the dissipative particle is determined as 

F = T In (2 sinh ^) - T J] In (2 sinh (12) 
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As a result, the free energy of the dissipative oscillator ( pi]) in the classic limit of high 
temperature can be expressed in the form |23[ 



The result (JT^) can be clarified in the following way. The mean squared displacement of the 
oscillator (|T^), according to the fluctuation-dissipation theorem, has the form 

<.^>= r^<^'>- (14) 



— oo 



where 



< ^ >- = 2 02^ • COth — 15 



The internal energy f/ is a sum of kinetic and potential energies 



1"°° duj f muj'^ mQ'^\ n 
U= ^{^ + ] <s^>u (16) 



—00 
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Eq. ^ follows also from the exact theory. The free energy F = U — TS, where the entropy is 
S = —dF/dT, can be found solving the differential equation what results in the expression 
(p!3|). In our over-damped case (^) one should put f2 = and in the limit of high temperature 
h'-y <t^T the free energy ( |13| ) is reduced to 

F=^ln^ (17) 
2 T ^ ^ 

When the viscous drag coefficient 7 tends to zero, the entropy of a free motion goes to 
infinity and the free energy decreases with the reduction of 7, according to Eq. O. As 



follows from Eqs. 11 and 10, the analogue of 7 in Eq. 17 is rjfii/ p and the free energy of 



fluid fluctuations is the sum of free energies of independent oscillators 

For fixed bodies the total number of fluid modes is not a fluctuating variable. The eigenvalues 
fii depend on a distance between bodies giving rise to an /^-dependence of F, what determines 
the hydrodynamic van der Waals interaction according to Eq. 0. 



IV. EIGENVALUES FOR THE PLANE GEOMETRY 



Below we find the van der Waals interaction energy for two infinite plates immersed into 
a fluid, parallel to {x, ?/}-plane and separated by the distance R. The 2;-axis is cut by planes 
at points z = ±R/2. The solution of Eq. ^with the definition = \J ^ — k"^ has the form 

72/ 



ayr, z) 



/d k ( ^ \ ^ 

■^^-p yAk cos zk^ + Bk sin zk^j exp {ikr) (19) 

where r = {x, y}. The solution of Eq. |^ can be presented as 

/<Pk 
- — -r (afc cosh kz + (3^ sinh kz) exp {ikr) (20) 

The boundary condition (|^) at 2; = ±i?/2 results in six equations for six independent 
variables a, (3, A^, B^, {kyA.j. — k^Ay), and {kyBj. — k^By). There are two sets of eigenvalues 
determined by the following relations 

fc, tanh Ktan = 0, K^tan h/ctanh = (21) 

2 2 2 2 ^ ' 



The solution of Eq. 21 can be written in the form 



2 2 

= qin + ^ttln , hz = q2n + ^a2n (22) 

Jri JrC 



where gi„ = (1 + 2n)7r/R and q2n = 2mi/R with n = 1,2, 3, ... 



V. FREE ENERGY OF THE FLUID BETWEEN TWO PLANES 

The eigenvalue ii = k"^ + k1 should be inserted into Eq. If the surface area of each 
plane is 5*0, the total free energy (|18]) can be expressed through the energy density per unit 
area of each plane Uvdw = F/Sq, which is given by the following relations 

U,dW = W, + W2 + W3, W,,2,3 = ^J ^/l,2,3 (23) 



In the expressions for 1^1,2,3 one has to omit i?- independent terms (surface energy) and 
i?-hnear terms (volume energy), which are not relevant for the interaction problem. In 
Eq. ^ the summation over n can be extended up to infinity. In the expression for /a the set 
{2n, {l + 2n)} is denoted as n and the diverging sum should be cut off at A; ~ nn/R ~ I/Rq. 
It is easy to see that at _R ^ _Ro Eqs. ^ and ^ produce the /^-dependent part of W2 

The quantity /s in Eq. ^ can be written as a sum = /j^"* + /g^"* , where 

n>l ^ ' n>l 

Comparing to Eq. in Eq. ^ the sum of constants is omitted since it leads to a i?-linear 
term according to the rule ^ R J dkz/2TT. The n-product in Eq. ^ is restricted by 

n < 7t/Rq where compression fluctuations become important. In reality, an effective cut off 

(2) 

of the contribution of transverse fluctuations in /g is smooth, what can be described by 
the relation 

ipix) is an even analytic function, ip{0) = 1, and ip{x) — s> at x — 00. Differentiation of 
Eq. ^ produces the expression 

dfP 1 ^ 1 <9 / 2 , fnRoW 1 



dR 2 ^ ^{nRo/R)R'^ dR \ \ R J J R 
The n-sum in Eq. ^ at i? ^ i?o can be substituted by an integration with the exponential 
accuracy df^"* /dR = 1/R — c/ Rq, where the numerical constant is 
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Finally, in the limit R^ Rq 



/f = -c^ + ln^ + Ci (31) 



where ci is a numerical constant. Eq. |3T] is derived neglecting the small parameter kRo. 
An account of finite kR^ leads only to i?-independent corrections of the result (|31|). Use of 
Eqs. ^7\, ^ and results in the following expression 



= I / " ^"^^^^^ ^^^^ 



After integration Eq. ^ reads 



VI. THE HYDRODYNAMIC VAN DER WAALS INTERACTION 



As follows from Eqs. and 24, 



The main contribution to Wi comes from k ^ 1/R , what allows to use the perturbation 
theory for ai„ = and = a(92n)- In the limit kR ^ 1 the function a{q) satisfies 

the equation 

g + — I tana = —k (35j 
R J 

The solution of Eq. |3^ can be presented in the form a{q) = ao + ai + 0^2 + ... where 

k 2kao{q) f , , ^1 ^ ^^ "?(^) 

ctQ = — arctan — , ai[q) = — ^ , Q;2(,q'j = I 1 + 7 arctan ■ ' 



(36) 

In the limit kR S> 1 in Eq. ^ one can change summation by integration according to the 
rule Yin ~* ^ I dq/T^-, what brings only an exponentially small error 
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The first (after a constant) non-vanisfiing term in /i is proportional to l/R^k"^. Tlie inte- 
gration in Eq. 37 yields 

57r2 



/i 



(38) 



The /c-integration in Eq. ^ is cut off at k <^ I/-R0 and results in the logarithmically big 
term 



Wi 



SttT , R 
m ■ 



(39) 



The total hydrodynamic van der Waals interaction energy per unit area of each plane with 



the logarithmic accuracy at Rq can be obtained from Eqs. P3l 26, and 



UvdW 



n T R 

"24 W 



(40) 



Eq. ^ is an exact analytic result for the case of two parallel infinite planes. The interaction 
(^) is mediated by by fluctuations with the wave length much shorter than a typical micron 
diameter of particles. In Eq. ^ the inter-plane distance R should exceed the length Rq, 
otherwise at smaller R longitudinal fluctuations also contribute to formation of interaction. 
The hydrodynamic van der Waals interaction is attractive due to suppression of dissipation 
by the restriction of motion of the fluctuating fluid which is provided by stick boundary 
conditions. This is equivalent to an effective reduction of viscous drag coefficients of fluid 
modes and hence, according to Eq. |1^, to a reduction of the interaction energy when the 
planes approaches each other. 

It is instructive to compare the obtained hydrodynamic van der Waals energy (|40|) with 
the analogous value m^^W '^iS^ frequency (visible light) electromagnetic fluctuations, 

what is the conventional van der Waals interaction With weakly dispersive in the 

visible spectrum dielectric constants of polystyrene walls 2.40 and water between them 1.77 
one can obtain 



he 



^-1.1x10-^ — 



(41) 
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The small numerical coefficient in Eq. ^ is an essential feature of the conventional (electro- 
magnetic) van der Waals interaction. This is exactly that what makes this interaction small 
at the micron distance. In contrast, the hydrodynamic van der Waals interaction ( ^0]) is not 
proportional to a small coefficient and at the distance R ^ 1/im 

^ ~ 10^ (42) 

„,em V / 

'^vdW 

-3, 



The estimate Rq ~ 10~ /im for water is used in Eq. Despite previous formulas refer 
to the case of infinite planes, the numerical coefficient in the hydrodynamic van der Waals 
interaction Uydw of two spheres also should not be expected to be small. As one can 
see in the case of conventional electromagnetic van der Waals potential, the same reasons 
are responsible for smallness of the numerical coefficients for both planar and double sphere 
geometries |2^. For micron size particles separated by micron distance, when a ~ i? ~ l/im. 



the interaction can be estimated as U^dw — UvdwR"^ resulting in the relation 

U^dw — -const X T (43) 
where the constant is of the order of unity. 

VII. THE PARADOX 

In case of conventional van der Waals interaction mediated by electromagnetic fluctu- 
ations the mean values of electric and magnetic fields are zero < E >=< H >= 0. The 
stress tensor (Th. for electromagnetic field is quadratic with respect to fields and hence the 
mean value < > is not zero. This makes the origin of the force due to electromagnetic 
fluctuations physically understandable. The situation with hydrodynamic fluctuation forces 
is different. The mean value of a solution of the linearized Navier-Stokes equation is zero 
< V >= (the same for a fluctuation part of p). The hydrodynamic stress tensor 

/ dvi dvk\ . 
a.. = ,^^ + ^J-p5.. (44) 
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is linear in fluctuation variables, its fluctuation part is zero, and a fluctuation force has 
to be zero in this approximation. A non-zero contribution to < > can result from 
the nonlinear terms in the Navier-Stokes equation neglected at the above approximation. 



Account of this non-linearity has been done in Ref. and a finite fluctuation force has 
been obtained. This result was shown to be incorrect in Ref. [^, where the exact mean 
value of the stress tensor was found to be zero on the base of exact non-linearity of 
the Navier-Stokes equations. A conclusion of Ref. was that a fluctuation interaction 
mediated by hydrodynamic fluctuations was impossible. This contrasts with the result (^0|). 
What is wrong? 

To understand the situation let us consider the linear chain of small particles, connected 
by elastic springs, shown in Fig. 1 and described by the dynamic equation 

= ^{Un+l + Un-l - 2Un) (45) 

where s is the sound velocity and ao is the period. Two spatially fixed big particles substitute 
small particles, as shown in Fig. 1. The system is elastic and the force acting on a big particle, 
placed on the site n is 

2 

T71S 

Fn = - Un-l) (46) 

% 

Here m is the mass of a a small particle. The energy of fluctuation motion of small particles 
Uvdw is determined by self frequencies in the system like in Eq. |12[ Two different positions of 
big particles, "natural" in Fig. la and "compressed" in Fig. lb, have identical self frequencies 
since the system is harmonic and hence U^^yy = ul^'Jyy. In this situation there is no van der 



Waals force, what can be understood since the mean value of the force (|46D, which is linear 
in displacement, is zero. On the other hand, for another "natural" position (Fig. Ic) self 
frequencies differ from those in Figs, la and lb and hence f/^^vi/ 7^ U^dw ^'^'^ ^^'^ "^^^ 
Waals force is non-zero. The reason for this is that a transition from the position (a) to the 
position (c) in Fig. 1 cannot occur within a harmonic approximation, one should destroy 
some harmonic springs and rearrange them again in a different way. The linear expression 
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for force (^61) is not valid to describe a transition from (a) to (c) and a real force is non-linear, 
what makes its average finite even for < Un >= 0. 

An analogous situation takes place in hydrodynamics. According to its derivation, Eq. ^ 
is valid only for discrete R = Nao, where ao is inter-atomic distance and is an integer 
number. The full dependence of Uy^w on R has a structure on the atomic scale corresponding 
to remove of discrete atomic layers from the inter-plane region. The hydrodynamic expres- 
sion for the stress tensor (|4^) is not valid at such short scale like Eq. ^ cannot describe 
breaking of harmonic bonds. In contrast to smooth van der Waals potential mediated by 
electromagnetic fluctuations, the interaction mediated by hydrodynamic ones has a struc- 
ture as a function of distance on the atomic scale superimposed on the smooth function 
(^). To some extend, this is analogous to observation of the structured interaction poten- 
tial [|26|, where the role of atoms was played by small particles. The resulting statement is 
that the hydrodynamic expression for the stress tensor (0) cannot be used for calculation 
of fluctuation forces since it becomes non-linear (and contributes to those forces) at short 
distances where the hydrodynamic approach is not valid. Hydrodynamic fluctuation forces 
should be calculated on base of energy like it was done above. The conclusion of Ref. on 
absence of hydrodynamic fluctuation forces based on use of the hydrodynamic stress tensor 
is incorrect. 



VIII. INTERACTION MEDIATED BY BODY FLUCTUATIONS 

Suppose two identical spheres of the radius a are separated by the center-to-center vector 
it and immersed into a fluid. There is no any macroscopic motion in the system except 
one caused by thermal fluctuations of body variables: linear velocities of spheres ui,2(t) 
and their angular velocities Qi 2{t). The hydrodynamic medium cannot fluctuate itself and 
their motion is only induced by fluctuations of body variables. The resulting free energy of 
fluctuations of body variables I{R) depends on the inter-body distance due to stick boundary 
conditions imposed on an induced fluid velocity. The energy I{R) is formed by fluctuations 
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of some characteristic frequency, specified below, for which fluctuations of the inter-body 
distance is relatively small. This allows to characterize the energy / by the inter-body 
distance R. In reality, fluctuation effects of a fluid (leading to the hydrodynamic van der 
Waals interaction Uvdw) and of body variables cannot be separated and they together result 
in the total fluctuation energy E{R) which is not the sum U^dw + I general case. Since 
the interaction is formed by by fluctuations of some typical frequency, this description is 
valid only for a time scale bigger than the inverse of that frequency. For shorter time scale 
the reduction of the interaction between bodies to a pairwise potential is not correct. 

There is a special case, when the total interaction energy of two particles, mediated both 
by fluid fluctuations v{f,t) and body fluctuations Ui^2{t), can be written as the sum 

E{R) = U^dw + I{R) (47) 

Here U^dw is contributed only by fluid fluctuations (bodies are at rest) and only body 
fluctuations are responsible for formation of / (fluid fluctuations are only induced). The 
separation of different mechanisms (^7]) is valid when the mass density po of particles is much 
bigger compared to the fluid density p. This limiting case of heavy particles is considered 
below what enables to understand the main physics and to draw a conclusion about the 
general case po ~ P- As follows from the linearized Navier-Stokes equation (||), a typical 
fluctuation frequency of the transverse fluid velocity (with divu = 0) is Uy = ujqo? / R^, where 

c^o = ^ (48) 

An estimate of a typical frequency of fluctuations of body variables follows from their equa- 
tion of motion. The linear velocities of two spheres obey the equations 

47r n dui 2 ^ , , s 

yPoa»-^ = fu (49) 

The forces exerted on spheres by the fluid are expressed by means of the friction tensors Co 
and Ci 

FU = -Co'<2 - Cl'ul, (50) 
14 



Eq. ^ with the estimate ~ 7r?7a gives a typical frequency due to the motion of body 
variables Uu = ujqp/pq. Under the condition 

^ « ^ (51) 

when the inequality p ^ Po holds, the typical frequency of body motion is much smaller 
than one for motion of the fluid u^. In the classic (only thermal fluctuations) case fluctuations 
of body variables, giving rise to I{R), are well separated in frequency from fluctuations of the 
fluid velocity resulting in Uydw{R)- On the short time scale cu"^ the distance R fluctuates 
slowly and can be considered as a constant. On the long time scale uj^^ the potential 
Uvdw{R) is already formed but, due to the condition d'^Uydw/dR^ <^ pQa^u"^, it violates 
weakly the dynamics of body fluctuations, which remains almost free. This leads to a 
possibility of separation of the two mechanisms in the total fluctuation energy (ff^). At 



frequencies u ^ tOu, according to [19|, one can omit the frequency dispersion of forces in 



Eq. ^ and the friction tensors, as shown in Ref. 1^^, have the forms at a ^ i? 



and 



For the motion described by Eqs. ^ ^ and |5^ the Reynolds number {pa/i])V < v? >, 
where < > is thermal mean-squared velocity, is small. The Ci-term in Eq. ^ results in 
a macroscopic force on the body (1) due to a macroscopic motion of the body (2). This is 



traditionally called the hydrodynamic interaction |2^ . 

The interaction I{R) is formed on the time scale l/cUu, where the typical frequency of 
fluctuations of body variables can be estimated as ~ 10~^sec for a micron size particle in 
water. On this time scale thermal fluctuations of the inter-body distance are much smaller 
than micron range separation distance. This justifies the concept of the pairwise potential 
I{R) on the time scale bigger than 1/uJu- Below we consider the inter-sphere distance R 
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satisfying inequalities 



p 



(54) 



which unify the condition (^TJ) of vahdity of Eq. ^ together with the vahdity condition of 
the forms ( |52D and (^) for friction tensors. Under the condition a <^ R one can ignore 
fluctuations of angular velocities. Let us introduce the two vectors Si and §2 proportional 
to the total and the relative displacement of bodies 

dS.2 



dt 



(55) 



As follows from Eqs. ^ and 0, the equations of motion for those displacements have 
the form (there is no summation on the index i) 



q2 ci f) Ci 

3 ^ ^1,2 , ^ ,i "^^1,2 



+ 6nriaX\ 2 



dt 







(56) 



where 



9a2 3a 



(57) 



The vector R is supposed to direct along the 2-axis. According to Eq. there are six 
independent fluctuating coordinates, which contribute to the free energy. Comparing Eq. ^ 
with the equation of motion (pTj), one can write down, analogously to Eq. the free energy 
of body fluctuations in the form 

where the interaction energy is 

3 

m = '^J2^n{X\Xl) (59) 

i=l 

Insert of Eq. |5^ into Eq. ^ produces the interaction energy due to body fluctuations 

HH)^%T^ (60) 
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The condition of applicability of Eq. ^is given by the criterion (^41). The interaction energy 
(|60|) grows up under reduction of the inter-body distance. To clarify the physical meaning 
of this result one can consider a motion of the body, when the second body is at rest for 
simplicity. The friction force exerted by the fluid on the moving body is enhanced by the 
resting body since it stops the on its surface. This increases the viscous drag coefficient of 
the moving body with the reduction of i?, what results, according to Eq. |l^ in the repulsive 
force caused by body fluctuations {dl/dR < 0). 



IX. DISCUSSION 

The total fluctuation energy of interaction of two particles in a fluid can be written as a 
sum of two different contributions (^) under the condition (|5TD when p ^ po- The repulsive 



part I{R) is found under the condition (|^). The van der Waals interaction Uydw{R) in 
Eq. ^ (attractive part) is not calculated exactly for the case of two spherical bodies, in this 
work only arguments for its estimation (^) are given. The flrst remarkable feature of found 
fluctuation interactions is that they exceed much the conventional van der Waals interaction 
mediated by electromagnetic fluctuations [|l^], which is approximately —(10^^ — 10^^)T at 



the surface separation of the order of l/im. The second feature of these interactions is 
their different characters: I{R) is repulsive and U^dwiR) is attractive. The third feature 
is essentially different behavior of these parts when two bodies are not in a bulk fluid but 
approach a wall |^ or are conflned between two walls 0. The interaction I{R) drops down 
while approaching a wall due to a reduction of the velocity fleld, transmitted from one body 
to another, what is caused by the stopping action of the wall (we do not consider in details 
how I{R) is modifled close to a wall). In contrast to this, the van der Waals part U^dwiR) 
decreases hardly near a wall since it is formed mainly by fluctuations of short wave lengths, 
of the order of Rq, which are weakly affected by approaching a wall. In the general case 
(po ~ p) the separation of contributions (0) does not hold but, like in the limit p <^ po, 
the contribution of body fluctuations provides a repulsive effect in the total interaction 
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energy E{R) and fluid fluctuations give the attractive van der Waals contribution. While 
approaching a wall the repulsive contribution reduces but the attractive is violated a little. 
This seems to be relevant for explanation of the observed appearance of attraction between 
two bodies when they approach a wall 0. In experiments 0, and the attraction 
potential was of the order of temperature what also correlates with the developed theory. 
For a detailed comparison with experiments further calculations based on the hydrodynamic 
approximation are required. In the work when two bodies approach a wall very close, 
the attraction disappears again. It can be explained by an increase of the Coulomb repulsion 
due to the electrostatic interaction through the dielectric wall instead of Debye's screening. 

The results of different experimental groups reveal a very stable feature of the particle 
interaction in a fluid: a value of the attractive potential is always of the order of temperature, 
what is a strong indication of a fluctuation nature of the attraction. In the works and 
particles were conflned between close glass walls, in the work only one glass wall was used, 
and in the works [0] and ^ particles were placed on the air-water interface. This variety of 
situations enables to reject some mechanism of attraction. For example, mechanisms related 
to a charge of the glass walls unlikely play a crucial role, since the air- water interface in the 
experiments ^ and essentially differs from the glass-water surface, nevertheless the value 
of the attraction potential was approximately the same like in the works [Q, [^], and [^. For 
particles on the air-water interface the concept of two types of interaction / and U^dw also 
works, but it needs some modiflcation. 

In summary, the new mechanism of interaction of particles in a fluid is proposed based 
on positional and rotational fluctuations of particles. This mechanism competes with the 
particle interaction mediated by hydrodynamic fluctuations of a surrounding fluid. These 
interactions are not small at a micron size inter-particle distance and their interplay cor- 
relates with the experimentally observed attraction between particles, which appears when 
they approach a wall inside a fluid. 

I would like to thank C.-K. Au, A. Ladd, M. Medina-Noyola, and R. Rajagopalan for 
valuable discussions. 
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FIGURES 

• • • Q • • • • • Q • • • (a) 



• • • Q** •••Q— • — • — • — (b) 



I O • • O • • • \ (c) 

FIG. 1. The linear chain of particles. Each horizontal segment between particles behaves like 
an elastic harmonic spring, (a) represents the "natural" positions of two attached particles (open 
circles), (b) is obtained from (a) only by a compression motion of attached particles without a 
destruction of harmonic bonds, (c) is another "natural" position obtained from (a) by destruction 
of harmonic bonds. 
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